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Graphs are closely related to quantum error-correcting codes: every stabilizer code is lo- 
cally equivalent to a graph code, and every codeword stabilized code can be described by 
a graph and a classical code. For the construction of good quantum codes of relatively 
large block length, concatenated quantum codes and their generalizations play an impor- 
tant role. We develop a systematic method for constructing concatenated quantum codes 
based on "graph concatenation" , where graphs representing the inner and outer codes are 
concatenated via a simple graph operation called "generalized local complementation." Our 
method applies to both binary and non-binary concatenated quantum codes as well as their 
generalizations. 



I. INTRODUCTION 

The discovery of quantum error-correcting codes (QECCs) and the theory of fault-tolerant 
quantum computation (FTQC) have greatly improved the long-term prospects for quantum com- 
munication and computation technology. This general QECC-FTQC framework leads to a re- 
markable threshold theorem, which indicates that noise likely poses no fundamental barrier to the 
performance of large-scale quantum computations 

Stabilizer codes, a quantum analogue of classical additive codes, are the most important class 
of QECCs 0,3. These codes have dominated the study of QECC-FTQC for the past ten years 
because of their simple construction based on Abelian groups. The recently introduced codeword 
stabilized (CWS) quantum codes framework 0-0] provides a unified way of constructing a larger 
class of quantum codes, both stabilizer and nonadditive codes. Based on the CWS framework, 
many nonadditive codes which outperform stabilizer codes in terms of coding parameters, have 
been constructed. 

Graphs are closely related to QECCs. It has been shown that every stabilizer code is local 
Clifford equivalent to a graph code 0, @] • The basic ingredients of a graph code are a graph and 
a finite Abelian group from which the code can explicitly be obtained [9(. Every CWS code also 
has a canonical form, where it can be fully characterized by a graph Q and a classical code C 0, 0]. 
So a CWS quantum code Q is usually denoted by Q = (G,C). When the classical code C is linear, 
Q is a graph code; therefore, graph codes, and hence stabilizer codes, are special cases of CWS 
codes. FIG. CD demonstrates the relationship between all quantum codes, CWS codes and graph 
(stabilizer) codes. 

For the construction of good QECCs of relatively large block length and good asymptotical 
performance, concatenated quantum codes and their generalizations play an important role 0, 0, 



id . 1 1 it ] . Combined with the CWS framework, families of good quantum codes, both stabilizer 



and nonadditive, have been constructed [l0|, UjJ] . Concatenated quantum codes also play a central 
role in FTQC, and the proof of the threshold theorem 0, 15|. Given the intimate relations 
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FIG. 1: Quantum codes 



between graphs and quantum codes, a question that arises naturally is whether there is a graphical 
description for concatenated quantum codes and their generalizations. Moreover, if there were such 
a description, for the case where both the inner and outer codes are CWS codes, the next question 
is whether the corresponding graph captures the "quantum nature" of the concatenated code. 

Previously, some related results on graph codes have been obtained. For instance, concatenation 
of graph codes may be described graphically by adding some auxiliary vertices. However, it remains 
unclear what the final graph after removing those auxiliary vertices will look like [16j |. The known 
examples of generalized concatenated codes only provide graphical descriptions in the case where 
the outer code is of a special form [rH ]. However, none of these previous works provides a 
general systematic graphical description for constructing concatenated quantum codes. Lack of 
such a description seems to indicate that using graphs to describe quantum codes was a very 
restricted approach. This issue will be addressed in the present work by developing a systematic 
method for constructing concatenated quantum codes based on a graph operation called "graph 
concatenation." 

To be more precise, we construct the concatenated quantum code 

Qc = Q in C Qout, (1) 

where the inner code Q m = (Q\ n ,C- m ) and the outer code Q ou t = (£/out)Cout) are both CWS codes. 
We require C m to be linear, but C Q ut can be either linear or nonlinear. Since C m is linear, Q m is a 
graph (stabilizer) code. We can then denote the parameters of Q m by [[n, k, d]] p . For simplicity, 
throughout the paper we assume that p is a prime number. When Q m encodes k qupits, the 
corresponding outer code Q ou t of length n' must be a subspace of the Hilbert space %® k n , i.e., we 
can denote the parameters of Q ou t by ((n',K',d'))k. When C out is linear, then Q out is a graph 
(stabilizer) code that can also be denoted by [[n' , k' , d']] p k, where K' = p kk> . 

We now state our main result. 
Main Result: The concatenated quantum code Q c can also be described as a CWS code, i.e., 

Qc = (Qc,Cc), (2) 
and Q c can be constructed via the following way: 

Qc = Qin E Qout 

— (Gin>Ci n ) C {Q outi C out) 

— {Gin C GoutiCin d Coutji (3) 

where C c = Ci n C C ou t is the classical concatenated code with the inner code Ci n and the outer code 
C ou t, and the graph concatenation Qi n C Qout i n Eq. (EJ) gives the graph Q c . And, we show that Q c 
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can be obtained by concatenating Qi n and Q ut via a simple graph operation called "generalized local 
complementation. " 

The main advantage of constructing concatenated quantum codes via Eq. ([3]) is that the "quan- 
tum part" of this construction is fully characterized by the graph concatenation Q m C Q OU f Pro- 
viding the rule for performing this graph concatenation, the problem of constructing concate- 
nated quantum codes becomes purely classical, i.e., constructing the classical concatenated code 
Ci n C C ou t- Despite the restriction that C; n must be linear, our method of graph concatenation can 
be applied to very general situations: both binary and non-binary concatenated quantum codes, 
and their generalizations. 

This paper is organized as follows. In Sec. HU we give a simple example which informally 
demonstrates the rule of graph concatenation via generalized local complementation. In Sec. IIIIl 
we review definitions of graph states, CWS codes, and graph codes. In Sec. \YV\ for a simple case 
that the inner code encodes only a single qupit (i.e., k = 1) and the outer code is also a graph code, 
we provide a description of graph concatenation based on the algebraic structure of stabilizers. We 
prove our main result in Sec. [V] In Sec. IVIl we discuss the application of our main result to the 
situation of the generalized concatenated quantum codes. A final discussion and conclusion is given 
in Sec. EH 



II. A SIMPLE EXAMPLE AND THE RULE 

In this section we give a simple example to demonstrate the idea of our main result given by 
Eq. ([3]). We first recall how to describe a CWS code by a graph and a classical code; then we 
demonstrate how to represent a concatenated quantum code as an encoding graph Q^ nc ^ with 
some auxiliary vertices. Finally, we show how to obtain the graph Q c of the concatenated code 
c ), via "generalized local complementation" and removal of the auxiliary vertices. 



A. The graph and the encoding circuit of a CWS code 

Let us start by taking the outer code to be a simple n = 3 binary CWS quantum code Q Q ut = 
({?out> Cout), wnere ^out is a triangle given in FIG. [2JA Q out defines a unique quantum stabilizer 
state, which we denote \ip)g out - C ou t is a classical binary code of length 3, and can be either linear 
or nonlinear. A basis of the CWS code Q Q ut can then be chosen as {Z Cout \ip)g out }, for all the 
codewords c out G C ou t 0]. 

If C out is linear, then Q out is a graph (stabilizer) code. The name "graph code" is chosen due to 
the fact that there is a graphical way to represent the code Q ut 5 which gives both the information 
of £out and C out 01- 

To show how to represent C ou t graphically, let us first recall the encoding circuit of Q ou t = 
({?out> Cout)- We use the standard quantum circuit notations, for instance as those given in For 
a CWS code, in general the encoding can be done by first performing a classical encoder C^f^ 
which encodes the classical code C ou t and then a graph encoder G ou ™ c ^ which encodes the graph 
state corresponding to the graph Q out 0| as shown in the top left circuit of FIG. El Here qi,q2,Qs 
denote qubits 1,2,3 in FIG.[2A, respectively. 

Let us now take C ou t = {000, 111} which is linear and gives a [[3, 1, 1]] stabilizer code. In this 
case, the classical encoder C^°^ which encodes 

000, 1 111 (4) 
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can be implemented by adding an input qubit go and performing controlled-NOT gates with control 
qubit go an d target qubits gi , g2 , <Z3 , followed by measuring go i n the Pauli X basis (which can be 
done by applying a Hadamard gate on go and then measuring in the Pauli Z basis) as shown in 
the top right circuit of FIG. El Throughout the paper we always assume that we get the desired 
measurement outcome (if not, we just need to perform some local Pauli operations according to 
the actual measurement outcome). 



ki>- 








M — 


s~i{enc} 
^out 




s~i{enc} 
^ out 


\Qs) — 









m 

|9i>-& 
1 92 

1 93 



H 



-e- 



G 



{enc} 
out 



9o 



|9i>"^ 

1 92) — 

1 93) — 



90 



H 



H 



-0- 




FIG. 3: Encoding circuit for the [[3, 1, 1]] outer code. 



The graph encoder G^f^ consists of three Hadamard gates on gi , g2 , g3 and three controlled- 
Z gates between them (controlled-Z gates are applied whenever the corresponding vertices are 
adjacent in graph FIG. 0A.), as shown in the bottom left circuit of FIG. [3j Now it is clear that we 
can "move" the classical encoder C^°^ to the right of the graph encoder G^ c ^ by replacing each 
controlled-NOT by a controlled-Z, as shown in the bottom right circuit of FIG. El 

In the following we use the convention to modify the encoding circuit by applying a Hadamard 
gate on the auxiliary qubit go before applying the classical encoder as shown by the left circuit of 
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FIG. \E\ This modification can be viewed as "a basis change" of the input qubit qo, i.e., what the 
"classical encoder" C^ c ^ does is then 

+ -> 000, - -> 111, (5) 
where ± are the labels of the quantum states 

|±> = -L(|0)±|1)). (6) 

This change of basis yields a non-classical encoding circuit, yet we know that it does not make 
any difference for the quantum code because by this new encoding circuit we obtain the same code 
space as before. We adopt this convention throughout the paper: for any CWS code Q = (G,C), we 
always assume that the "classical encoder" C^ enc ^ maps "classical strings" in the {|+), |— )} basis 
to "classical codewords" in the {|0),|1)} basis. We will see later that this convention naturally 
leads to a simple rule for graph concatenation. 

Moreover, for any [[n,k, d]] CWS code Q = (G,C) with linear C (i.e., Q is a graph code), the 
encoding of Q can be applied by first performing the graph encoder G^ enc ^ , and then the classical 
encoder C^ enc ^ as follows: use k input qubits; apply Hadamard on each of the k qubits; replace 
each controlled-NOT gate performed in the original classical encoder C^ enc ^ with a controlled-Z 
gate; finally measure each of the k auxiliary qubits in the Pauli X basis. 

This encoding circuit can be represented graphically: add the k input qubits as k new vertices 
to the graph G; whenever a controlled- Z is applied in the encoding circuit between an input vertex 
v and a vertex v' of Q, add an edge between them [||. The corresponding graph representing the 
graph code Q = {Q,C) is denoted by Q c . 

Therefore, for the outer code Q out = (G OVL t,C ou t) with graph Q out given in FIG. [2]A., where 
C ou t = {000, 111} is linear, we can insert the input qubit qo as a new vertex (denoted by 0) to 
the graph Q (the middle white vertex in FIG. 03). We then add the edges between and 1,2,3 
according to the encoding circuit given by the bottom right circuit in FIG. [3] (see FIG. |2j3) . This 
graph is then denoted by Q on X ■ There are two types of vertices in Gq^ 1 '■ the input vertices (the 
middle white vertex); and the output vertices (vertices 1,2,3). 



B. The encoding graph of a concatenated quantum code 

Now we consider the inner code Q m = (Gi n ,Ci n ). Notice that due to our restriction for Eq. ([3]), 
Ci n must be linear. So Qj n is a graph code and has a graph representation <7 ir j n . For simplicity we 
take Q; n be a [[2, 1, 1]] stabilizer code, which is represented by the graph of FIG.[4]A on the vertices 
1,4,5. (The subgraph of the vertices 4,5 represents Qi Q , and 1 is the input qubit describing the 
classical encoder of C; n ; hence, C m = {00, 11}.) 

To construct the concatenated code Q c = <2i n C Qout, since the outer code has length n' = 3, 
we must take three copies of <5; n , to encode qubits 1,2,3 as shown in FIG. 0]A.. The graphical 
representation of the concatenation procedure is shown in FIG. Hf3- Here, in the outer code, the 
middle white vertex is encoded into vertices 1,2,3. Then each of these vertices is encoded using 
the inner code: vertex 1 into vertices 4, 5; vertex 2 into vertices 6, 7; and vertex 3 into vertices 8, 9. 

We call this graphical representation of the concatenated code Q c with a linear C out the encoding 
graph of Q c and denote it by g^tienc} have three types of vertices in g c ^{ enc } : fa e i n p U t 
vertices (the middle white vertex in our example); auxiliary vertices which are in the subgraph 
Gout (vertices 1, 2, 3); and output vertices which are in the copies of Q m (vertices 4, 5, 6, 7, 8, 9). In 
general, if C ou t is nonlinear, similarly we can have an encoding graph of Q c and denote it by Gq HC \ 



FIG. 4: Concatenated graph code. 



which in our example is the subgraph of FIG. HJ3 without the middle white vertex. Therefore, we 
only have two types of vertices in Gq^'- the auxiliary vertices (vertices 1,2,3); and the output 
vertices (vertices 4, 5, 6, 7, 8, 9). 

The encoding circuit of the concatenated code Q c is given by the right circuit in FIG. where 
G^ nc ^ denotes the graph encoder for the graph of the inner code Q m . To obtain this encoding 
circuit, we should recall our convention of adding a Hadamard gate before performing the classical 
encoder. 
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FIG. 5: Encoding circuit for the concatenated code with linear outer code. 



In general, if C ou t is nonlinear, the encoding circuit of the concatenated code Q c is given by the 
right circuit in FIG.[H where G^ nc ^ denotes the graph encoder for the graph of the inner code Q m . 
Again, note that we add a Hadamard gate before performing the classical encoder. Also, we need 
to keep in mind that the "classical encoder" C^"^ maps "classical strings" in the {|+), |— )} basis 
to "classical codewords" in the {|0), |1)} basis. 
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FIG. 6: Encoding circuit for the concatenated code with a general outer code. 



C. The rule of the generalized local complementation for graph concatenation 



As shown in Sec. IIIB( given a concatenated quantum code Q c = Q m C Qout with a CWS 
outer code Q out = (</out,Cout) and a g ra P n inner code Qout = (^out,C out ), it is easy to get the 
encoding graph Qq"^. We claim (and will show later in Sec. |Vj that the concatenated quantum 
code Q c = Qi n C Qout can also be described as a CWS code Q c = (G C ,C C ). Therefore, the real 
description that we want for the concatenated code Q c is a graph Q c and a classical code C c such 
that Q c = (Q C ,C C ). Also, we want the classical code be given by the "classical concatenation" of 
the classical code of the inner and outer code, i.e., C c = C m C C ou t5 so that the quantum part 
can be fully taken care of by the graph concatenation Q c = Q in C Gout- Furthermore, this graph 
concatenation should be given by some simple graph operations on the encoding graph Qq* 10 ^ , i.e., 
we want a general rule which gives 

4? c} -> So, (7) 



or 

gC^{enc}^ g e % (g) 

if the outer code is also a graph code. As discussed in the main result, such a general rule does 
exist and we call it "generalized local complementation." 

In this section, we demonstrate the rule of generalized local complementation for graph con- 
catenation by a simple example, starting from the encoding graph given by FIG.|3£>. Keep in mind 
that we want 



C c = C in c C out = {00,11} C {000,111} = {00 0000,111111}. (9) 

Remark 1 To obtain the graph Q c (or Gc c ) from FIG. \J$3, a naive way is to calculate the stabilizer 
state of the output vertices after performing Pauli X measurements on all the input and the 
auxiliary vertices in the encoding circuit (given by the right graph of FIG. or the right graph 
of FIG. [5j), and then to represent it as a graph state (or a graph code). Notice that in general it 
might not be possible to represent the very code as a graph Q c (or Qc c ) does not necessarily exist. 
Indeed, any stabilizer state is local Clifford equivalent to a graph state (which is not necessarily 
unique), any any stabilizer code is local Clifford equivalent to a graph code, and any CWS code is 
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local Clifford equivalent to a standard form given by a graph and a classical code. However, for a 
general CWS code, those local Clifford operations transform both the graph and the classical code 
J^j/. Therefore, it is not clear that such a graph Q c exists such that the classical code is obtained by 
concatenation, i.e., C c = Ci n C C ou t- 

Now we specify the rule of generalized local complementation: given a graph Q, for any vertex i, 
denote the set of its adjacent vertices by N(i). Also let S be a subset of vertices disjoint from N(i). 
A generalized local complementation on i with respect to S is to replace the bipartite subgraph 
induced on N(i) U S with its complement. 

For an example, the generalized local complementation of the graph shown by FIG. on 
vertex 1 with respect to S = {6, 7} results in the graph shown by FIG. [7£>. Here N(l) = {2, 3, 4}. 
The generalized local complementation replaces the bipartite subgraph of vertices {2, 3, 4, 6, 7} and 
edges {(3, 6), (4, 7)} with its complement (i.e. another bipartite subgraph of vertices {2,3,4,6,7} 
and edges {(2, 6), (4, 6), (2, 7), (3, 7)}). 
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FIG. 7: Generalized local complementation 

Now we are ready to specify the rule of graph concatenation in terms of generalized local 
complementation. (Recall that our goal is to obtain Q c from Qq , or Q^ c from g^ ut ^ enc ' 

Procedure 1 (Graph Concatenation via Generalized Local Complementation) 

1. Given the graph Qq C ^ (or g < £° nt { en °}) ) j or ea ch auxiliary vertex i, define Si to be the set of 
all output vertices which are adjacent to i. 

2. For each auxiliary vertex i, delete all the edges which connect i to vertices in Si- 

3. For each auxiliary vertex i, perform generalized local complementation on i with respect to 
Si. Note that the order in which we apply those generalized local complementations does not 
matter since the whole procedure on all auxiliary vertices finally gives the same graph. 

4- Remove all the auxiliary vertices. 

To demonstrate the above rules, let us start from the encoding graph given in FIG. HJ3 for the 
concatenated quantum code with the outer code [[3, 1, 1]] given in FIG. [2)3 and the inner code 
[[2, 1, 1]] given in FIG. HA. For convenience we redraw FIG. E£> in FIG. [BJA.. Now from FIG. [8]A. we 
get Si = {4,5}, 52 = {6,7}, and S3 = {8,9}. Deleting all the edges which connect each auxiliary 
vertex i and vertices in Si (for i = 1,2,3) results in FIG. [8)3. Performing local complementation 
on auxiliary vertex 1 with respect to Si = {4, 5} leads to FIG. [SC, where we use dashed blue lines 
to show the edges that we add between output and auxiliary vertices, and solid black lines to show 
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FIG. 8: Generalized local complementation for graph concatenation 

other edges. Performing local complementation in FIG. [5TJ on auxiliary vertex 2 with respect to 
S2 = {6, 7} leads to FIG. [Hp, and performing local complementation in FIG. [Hp on auxiliary vertex 
3 with respect to £3 = {8, 9} leads to FIG. [Hp. Removing all the auxiliary vertices in FIG. [HE gives 
FIG. EF, which is the graph Cfc . 

From FIG. [Hp one can easily see that the rule for concatenation of the classical codes given in 
Eq. © holds. In general, the outer classical code C out is nonlinear, so we do not have the input 
vertices in the encoding graph of the concatenated code. However, we can still go through the whole 
procedure of the generalized local complementations on auxiliary vertices to obtain the graph Q c . 
In our example this procedure is demonstrated by subgraphs of FIG. [8j\ through FIG. [Hp without 
the middle white vertex. 

III. GRAPH STATES, CWS CODES, AND GRAPH CODES 

In this section we review the stabilizer formalism to fix our notation especially in the non-binary 
case; then we define CWS codes, graph codes, and finally describe their encoding circuits. 

A. The generalized Pauli group 

Let p be a prime number and ¥ p be the field of p elements. A qupit is a p-level quantum system 
whose Hilbert space is represented by the orthonormal basis {|r) : r £ F p } = {|0), |1), . . . , \p— 1)}. 
Let uj = e 27rl//p be a p-th root of unity. The (generalized) Pauli matrices X and Z are defined as 
follows. 



X\r) = \r + 1 mod p), 



(10) 
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Z\r) 



uj r) 



(11) 



It is clear that X p = Z p = I, and then we can consider the operators X a and Z where a, b E ¥ p . 
We have Z b X a = uj ab X a Z b ; therefore, X a Z b and X a ' Z y commute iff ab' - ba' = (see e.g. Q 
for more details). 

The group generated by the (generalized) Pauli matrices X and Z is {uo c X a Z b : a, b, c G F p } 
and is called the (generalized) Pauli group. Notice that, in the binary case (p = 2) the Pauli group 
is generated by Pauli matrices a x and a z together with il {i = \f— T). 

Let n be an arbitrary positive integer. For vectors a = (ai, . . . , a n ) and b = (6i, . . . , 6 ra ) in F™ 
define 



(12) 



and 



Z b = Z bl 



Z 



b„ 



(13) 



Again, two Pauli matrices X a Z h and X a ' Z h ' commute if and only if ab' — a'b = 0, where cd = 
c\d\ + ■ ■ ■ + c n d n is the usual inner product on F™. 

For simplicity, a Pauli operator X a Z h is denoted by the vector (a | b) of length 2n. Thus two 
Pauli operators X a Z h and X a Z h commute iff their corresponding vectors are orthogonal with 
respect to the "symplectic inner product" defined by 



(a | b) * (a' | b') = ab' - a'b. 



(14) 



B. Stabilizer states 

It is easy to see that for a Pauli matrix g = co c X a Z h , g p = I. (In the case p = 2, the statement 
might only be true after replacing G by ig in order to get g 2 = I. Having this in mind, there is no 
true difference between the binary and non-binary case in the rest of the paper.) Therefore, the 
eigenvalues of g are all p-th root of unity. In fact, if (a | b) is non-zero, then for any i, to 1 is an 
eigenvalue of g, and the multiplicity of each of these p eigenvalues is equal to p n ~ l [20j]. 

1 i 1 k i k 

Now suppose gi = uj Cl X a Z , . . . , = uj Ck X a Z are k Pauli matrices which pairwise com- 
mute and such that the subgroup generated by any k — 1 of them does not contain the other 
one. Additionally, we require that the group generated by g%, . . . , g^ does not contain a non-trivial 
multiple of identity. Since gi,...,gk commute, they can be diagonalized simultaneously. 

Lemma 1 The common eigenspace of all gi's with eigenvalue 1 is a p n ~ k -dimensional subspace. 



This lemma is a well-known fact in the binary case [l(, and a proof for the non-binary case can 
be found in 20]. 

Representing the operators g\, . . . , gk by the vectors of length 2n, we obtain the k x (2n) matrix 



M 



b 1 \ 
b 2 



(15) 



\ a 1 * 



b k / 
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of rank k (because gt is not in the subgroup generated by the rest of g/s). The rows of M are 
mutually orthogonal with respect to the symplectic inner product (see Eq. (114D ). 

If we consider n generators, or equivalently an n x (2n) full-rank self-orthogonal matrix M, 
Lemma Q] implies that the common eigenspace of all gi 's with eigenvalue 1 is a one-dimensional 
subspace. Hence there is a unique (up to a scaler) non-zero vector \ip) such that gi\ip) = \ip). In 
fact, if we consider the group S generated by the g^s, for any h 6 S we have h\ip) = \ip). The 
group S, which is a maximal Abelian subgroup of the Pauli group modulo its center, is called a 
stabilizer group, and the state \tp) is called a stabilizer state. 

Notice that for a stabilizer state \ip), its stabilizer group S is unique; however, {gi, . . . ,g n } is 
just some set of generators of S. Suppose {hi, . . . , h n } is another generating set of S. Then for 
any i there is Uij £ ¥ p such that h{ = g^ 11 ■ ■ ■ g^ in . As a result, the vector corresponding to hi is 
equal to (un, . . . , u- m )M. 

Lemma 2 Any set of generators of the stabilizer group S with k generators can be represented by 
a matrix UM, where U is an invertible k x k matrix. 



C. Clifford group 

The Clifford group is the normalizer of the Pauli group. In the binary case, it is well-known that 
the Clifford group is generated by the Hadamard gate, the phase gate, and the controlled-NOT 
gate (l7j ]. A characterization of the Clifford group in the non-binary case can be found in (l8| . 

Clifford operators are important in the stabilizer formalism because they send any stabilizer 
state to a stabilizer state. Suppose is a stabilizer state with the stabilizer group S. Also, let 
L be a Clifford operator. For any g S S we have (LgL^)L\ip) = L\ip). On the other hand, LgL^ 
is in LSL^ which is a subgroup of the Pauli group since L is a Clifford operator. In fact, LSL^ 
is a maximal Abelian subgroup of the Pauli group whose corresponding stabilizer state is L\ip). 
Therefore, Clifford operators send stabilizer states to stabilizer states. 



Based on the characterization of the Clifford group [171. Il8l|. for any two stabilizer states 
and there is a Clifford operator L such that L\ip) = However, it does not mean that all 
the stabilizer states are the same in the point of view of quantum coding theory since the operator 
L may completely change the entanglement of a state. But if we assume that L = Li • • • ® L n 
is a local Clifford operator (L is the tensor product of n one-qupit Clifford operators), then the 
entanglement of and L\ip) are the same. Based on this idea, two stabilizer states are called 
"local Clifford equivalent" if they are equivalent under the action of the local Clifford group. 

For the encoding circuits, we need only two Clifford operators which we describe next. Define 
the vector 

|f> = -^|>-"|->, ( 16 ) 

V P 3=0 

for any r £ ¥ p . \r) is an eigenvector of X, i.e., X\r) = uj r \r), and {|0), . . . , \p — 1)} is an orthonormal 
basis. Therefore, the operator 

H\r) = \r), (17) 

which is called the (generalized) Hadamard gate, is unitary. By definition HXH^ = Z. Also, it is 
easy to see that HZH^ = Jft. Hence, H is in the Clifford group. Using the above relations the 
proof of the following lemma is easy. 
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Lemma 3 Suppose \tp) is a stabilizer state whose stabilizer group is represented by the n x (2n) 
matrix M. Thus, the matrix representation of the stabilizer state Hi\tp) (Hadamard gate is applied 
on the i-th qupit) is obtained from M by exchanging the i-th and (n + i)-th columns and then 
multiplying the i-th column by —1. 

The next operator is a two-qupit gate which is called controlled- Z and is defined by 

C z \r)\s) = \r)Z r \s) = Z s \r)\s) = u rs \r)\s). (18) 



We have 

C Z X ® IC\ 
C Z I ® XC\ 

c z z® ic\ 
C z l® zc\ 

and thus by definition C z is in the Clifford group. 



= X®Z, 

= z®x, 

= Z®I, 

= I®Z, (19) 



Lemma 4 Suppose \ip) is a stabilizer state whose stabilizer group is represented by the n x (2n) 
matrix M. Thus, the matrix representation of the stabilizer state C' z \ip) (the controlled-Z gate is 
applied on the i-th and j-th qupits) is obtained from M by adding column i to column n + j, and 
column j to column n + i. 



D. Graph states 

In the following we consider graphs whose edges are labeled by non-zero elements of ¥ p . Con- 
sidering the adjacency matrix of a graph Q, we can represent it by a symmetric matrix over ¥ p with 
zero diagonal. Suppose G is such a matrix of size n x n. Then M = (J n | G) is a full-rank n x (2n) 
matrix, and all of its rows are mutually orthogonal with respect to the symplectic inner product; 
therefore, M represents a stabilizer group which corresponds to a stabilizer state. Such a stabilizer 
state is called a graph state, which we denote by \ip)g- It is well-known that any stabilizer state 
is local Clifford equivalent to a graph state so to study the properties of stabilizer states it is 
sufficient to restrict ourselves to graph states. 

Graph states can be generated easily using only Hadamard and controlled-Z gates, 

Lemma 5 The graph state corresponding to the graph with adjacency matrix G = (gij) on n 
vertices can be generated by the following circuit. Prepare n qupits in the state |0), apply H* to 
every one of them, and then for any i,j apply C z tj on qupits i and j. 

Proof: The initial state of the n qupits is |0) • • • |0), which is a stabilizer state with the stabilizer 
group {Z a : a £ F™}. This stabilizer group corresponds to the matrix Mq = (0 n | I n ). According 
to Lemma [3l after applying W gates the matrix Mq will be changed to M\ = (I n \ 0). Also, by 
Lemma IH applying C z lJ on qupits i and j corresponds to adding columns i and j multiplied by 
gij to columns n + j and n + i, respectively. Since the first block of Mi is identity, this operation 
is the same as to add g^ to the entries ij and ji of the second block. Therefore, at the end we 
obtain the matrix M2 = (I n \G). □ 
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E. Measurement on graph states 



Suppose we have a graph state \ip)g which corresponds to the graph Q with adjacency matrix 
G, and we measure its (say) last qupit in the standard basis and get |0). We claim that the state 
after the measurement (without the measured qubit) is also a graph state whose corresponding 
graph is obtained from Q by removing the last vertex. To see this fact precisely notice that since 
(l n | G) represents the stabilizer group of \ip)g, for any i, we have X ei Z Sl \ip)g = \ifi)g, where all 
coordinates of ej are except the i-th one which is 1, and gj is the i-th row of G. Let 

p-1 

lV>)s = 2>r#r)|r), (20) 

r=0 

and for 1 < i < n — 1 let and be the vectors of length n — 1 obtained from gj and e[, 
respectively, by deleting the last coordinate. Thus we have 

p— l p— l 

\^) g = X^Z^\i)) g = ^a r Z9™X<Z^\(t> r )\r) = a r (u r9 ™X<Z^\(j) r )^ \r). (21) 

r=0 r=Q 

As a result X e iZ s 'i\<f)o) = \4>o), which means that \4>o) is a stabilizer state with the stabilizer group 
generated by X e >Z s % 1 < i < n — 1, and the matrix representation of these generators is (I n -i \ G') 
where G' is the adjacency matrix of the graph obtained from Q by removing its last vertex. 



F. CWS codes and graph codes 

A CWS code ((n, K,d)) p is described by a graph Q with n vertices and edges labeled by F p , 
together with a classical code C which consists of K vectors in Such a code is denoted by 
Q. (G.C) 00. 

If the classical code C is linear, then Q is a graph (stabilizer) code [j4a]- The parameters of 
such a graph code Q = (Q, C) are [[n, k, d]] p , where the classical code C consists of K = p k vectors 
in ¥p that are indexed by the elements of This [[n, k, d]] p graph code encodes k qupits into n 
qupits in the following way. Suppose \ip)g is the graph state corresponding to Q. To encode a state 
of the form (g) • • • ® H^\r\ . . . r^) we first find the classical codeword a G C which is indexed by 
r\ . . . rfc, and then encode ® • • • ® H^\r\ . . . r^) into Z a \ijj)g 2l|. Since C is a linear code, it is 



a linear subspace of F" We can then represent C by k basis vectors otx, . . . , ct^. In this case, the 
state Ft ® • • • ® Eft\ri . . . r k ) is encoded into Z riai+ - +rkak \ip)g. 

The encoding circuit of a [[n, k, d]] graph code is simple, as shown in the following procedure. 

Procedure 2 (Encoding circuit for a graph code) 

1. First generate the graph state \ip)g using the circuit described in Lemma\^ 

2. For any 1 < i < k apply on qi, where qi, ■ ■ ■ ,qu are the qupits that we want to encode. 

3. For any 1 < j < n apply Cz lJ on qi and the j-th qupit of\ijj)g, where atj is the j-th coordinate 
ofoii. 

4. Apply H to gi,... ,q k - 



5. Measure q%, . . . ,q k in the computational basis. 
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For example, the encoding circuit of the graph code with a triangle graph and the classical code 
{000, 111} can be found in the left circuit of FIG. EJ (Notice that in the binary case = H.) 

In general, for a graph code Q = (G,C) the encoding circuit can be represented graphically, 
and the corresponding graph is denoted by Q . consider the graph Q, for any 1 < i < k add a 
vertex (input vertices), attach it to the vertices of Q (called the output vertices), and label the 
edge between this vertex and the j-th vertex of Q by ctij. For example, FIG. [2fB gives the graph 
Q c , where Q is a triangle and C = {000, 111}. 

Remark 2 The encoding circuit corresponding to the graph code with graphical representation Q c 
is related to the circuit that generates the graph state \ip)gc corresponding to the graph Q c . To see 
this, notice that the steps 1,2, 3 in Procedure® indeed give such a graph encoder. 



To find the logical X and Z operators of an additive graph code we first describe the stabilizer 
group of the logical |0 . . . 0)l state. Notice that 



|0...0) 



1 



■■ E 

ri,—,r k 



H^\ n ...r k ), 



(22) 



and then 



|0...0) L 



(23) 



Therefore, all operators Z ai are in the stabilizer group of |0...0)x, 3 and the logical Z operators 
are described by the rows of the matrix 



(24) 



Since the vectors ctj are linearly independent, without loss of generality (by a change of basis for 
the classical code and reordering the qupits), we may assume that the first block of the second part 
of this matrix is Ik- So we assume that the matrix [Ik A), where A is of size k x (n — k) describes 
a basis for C, and the logical Z operators are 



( 


"1 \ 


V 


c*fc / 



(0 I I k A) 



Assume that 



G 



d B 
B T G 2 



(25) 



(26) 



where G\, G2, and B are of size k x k, (n — k) x (n — k) , and k x (n — k), respectively. Then the 
stabilizer group of the state \ip)g is represented by 



h 

I n -k 



Gi B 

B T G, 



(27) 



Now note that for any 1 < i,j < k, (Z ai )(X e i Z s j) = uj s » ' {X e i Z s i){Z ai ), where 5ij is the Kronecker 
delta function. On the other hand, the code space is invariant under X e i Z s i . Therefore, the logical 
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X operators can be described by the matrix 

(I k 0|d B), 

Also, it is not hard to see that the Pauli matrices corresponding to the rows of 

{-A T I n . k | -A T Gi + B T -A T B + G 2 ), 



(28) 



(29) 



commute with both logical X and logical Z operators. Therefore, the additive graph code Q is 
described by the stabilizer group 





~A^ In—k 



h A 
-A T G 1 + B T -A T B + G 2 



logical Z operators 



and logical X operators 



(30) 



(31) 



(32) 



Z=(00\ I k A), 

X = {l h | Gi B). 
G. Summary of notations 



Before going into the detailed proof of the main result, we summarize our notation. Let Q = 
(G,C) be a CWS code. If C is linear, then Q is a graph code, where the code has a graphical 
representation denoted by Q c . The concatenation of two CWS quantum codes Qj n = (<5m,Cm) and 
Qout = (£out,C ou t) is denoted by Q c = Q in C Q ou t- 

See Table I for the rest of notations. 



c 


the classical code 


c 


the generator matrix of the classical code C, if C is linear 


(j{enc} 


the encoder of the classical code C 


Q 


the graph corresponding to the graph state \ip)g 


G 


the adjacency matrix of the graph Q {G = (gij)) 


Q c 


the graph representing the graph code Q, if C is linear 


Q{enc} 


the encoding circuit of the graph Q 


p{enc} 

y Qc 


the encoding graph of the concatenated code Q c 


f-,Cout{enc} 

y Qc 


the encoding graph of the concatenated code Q c , if C ou t is linear 



TABLE I: Notations 



Most of the notations have already been given in Sec. IIIIl except for and Qq 

are discussed in Sec. Ill Bl and will be explained in more details in Sec. |V] 



Cout {enc} 



which 



16 



IV. CONCATENATION OF GRAPH CODES 



In this section, we prove our main result in a simple case, where the inner code encodes only a 
single qupit and the outer code is a graph code. In this situation, we can algebraically obtain the 
graph and the classical code of the concatenated code using the stabilizer formalism. Although we 
will prove our main result in the general case in Sec. [V] we believe that the proof given in this 
section is easily accessible to those who are familiar with the stabilizer formalism. 

Suppose the inner code Q m = {Gw,C\ n ) encodes only a single qupit, i.e., Q m is an [[n, l,d]] p 
code. Then from the discussion in Sec. IIIIF[ it follows that 



G in = ( y ° T y H , ) , (33) 



and 

C in = {0,(l b)}, (34) 

i.e., 

Cin = (1 b), (35) 

where both y and b are vectors of length n — 1. (Notice that Q m encodes one qupit; thus, C- m is 
one-dimensional. ) 

Since Q m encodes one qupit, the corresponding outer code Q out = (<5 ut, Cout) is an ((n', K\ d')) p 
code. In this section we assume that C out is linear, so Q Q ut is a graph code with parameters 
[[n', k', d']] p , where K' = p k . Then from the discussion in Sec. IIIIFl we have 

(36) 

and the rows of 

C out = {Iy A) (37) 

form a basis for C out . 

Thus by Eqs. ([30]) -(|32|) the stabilizer group of Q m is 



its logical operator Z is given by 





In-l 



1 b 

y T - h T y + H > 



(38) 



Z- m = (0 I 1 b) , (39) 

and its logical operator X is given by 

X\n = (l | y) ■ (40) 
The concatenated code Q c = Q m □ Q ou t is a quantum code which encodes k' qupits into nn' 
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qupits as follows; it first encodes kl qupits into n' qupits using Q ut> and then encodes any of the 
n' qupits into n qupits based on Q in . 

The main result of this section is given by the following theorem, which states that the con- 
catenated code Q c = Qi n C Qout is also a graph code. The corresponding adjacency matrix of the 
graph and the generator matrix of the classical code can be computed directly from the adjacency 
matrices and the generator matrices of the inner and outer codes. 

Theorem 1 Suppose Q out = (Gout,C ou t) and Qin = {Qin,C in ) are [[n', k', d']] p and [[n,k,d]] p graph 
codes, respectively, (where k = 1) as described by Eqs. ^SSfy-fflfy. Then the concatenated code 
Qc = Qin C Qout = (Ga Cc) is a graph code described by the graph Q c with adjacency matrix 



G c = G m ® I n , + I (1 b) ® G ouU (41) 




and the classical code with generator matrix 

C c = (l b)®(J fc / A), (42) 
i.e., the classical code is obtained by concatenation as well: 

Cc — C{ n rz C ut- (43) 



Proof: Let us first show that a basis of C c is described by Eq. ([42]). To find Z c , the logical Z 
operators of Q c , we should first consider the logical Z operators of Q ut> and then replace any Pauli 
matrix X and Z of those operators with the logical X and Z operators of Q; n . For example, if the 
logical Z operator acting on the first encoded qupit in Q out is Z^ 1 ' 1 ' '—' ', the logical Z operator 
acting on the first qupit in Q c is Z'(( 1,b ) , ( 1 > b )> 0- " ) because by Eq. (j39|) the logical Z operator of 
Qi n is Z^' h \ Therefore, by changing the order of qubits we can represent this Pauli matrix by 
the vector ^0 | (1 b) ® (1, 1,0,..., 0)j , where the zero before the vertical line is actually a zero 
vector. Now since the logical Z operators of Q Q ut are represented by rows of Eq. (J3JJ) , we have 



Z c =(0\ (1 b) ® {Iy A) ) . (44) 

Equivalently, (1 b) ® (J fc / A) is a basis for the linear code C c . 
Analogously, we compute for the logical X operators of Q c : 

X c = ((1 0) 8) {I k , 0) | (0 y) {I W 0) + (1 b) (8) {G x B)) (45) 
4 I Gi B b®(G!5)+y® {I k , 0)) . (46) 



It remains to compute the stabilizer group. The first n 1 rows of S c are obtained from rows of 
S ut (Eq. (|30p ) by replacing any X and Z with the logical X and Z of the inner code. For the next 
(n — l)n' rows note that, g<i, . . . ,g n which are the Pauli matrices corresponding to the rows 2, . . . , n 
of e>j n commute with the logical X and Z of the inner code. In fact, they are in the stabilizer group 
of the code space (spanned by the states |0), . . . , \p — 1) states). Now since we replace any of the 
n' qupits of Q Q ut with a state in the code space of Q m , each block of n qupits in Q c should be 
stabilized by 52, • • • ,9n- As a result, S c = (M\N), where 
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M 



(1 0)(8 
V ("b T /, 



-A T /„,_ 





n'—k' 



n-1. 



and 



TV 



(0 y)<8> 







V 

+ (1 b) 









In'-k' 



o \ 










Jn'-k' 



(n-l)n' 



(47) 



4' ^4 



(48) 



h'y + H')®!^ 



Now to complete the proof of Theorem [T] it is sufficient to show that the stabilizer group, and 
the logical X and Z operators of the graph code described by Eqs. (jH]) and (jUJ) are given by 
Eqs. (I44f) - (f48|) . We compute these matrices using the construction given by Eqs. (f30l) - (|32j) . 

First of all, the classical part of the code is given by (If./ B), where B = (Ah® (Iy A)); 
therefore, the logical Z operators of the code are the same as Eq. (j44]). 

The block from of matrix G c of Eq. (|4ip is given by 



Ki W 
W T K 2 



(49) 



where K 1 = G 1 ,W=(B y ® {I k > 0) + b ® (G x B)) and 



K, 



Gi 





Jn'-k' 



+ b 3 



B 

G 2 



y ® (0 In'-fe') + b (8) {B 1 G 2 
H' ® I n > + b T b G 2 



(50) 



Hence, the logical X operator of the graph code is 

Vfc, I K x W) = (iy | K x B y® [I h , 0) + b ® {G x B)) 



(51) 



which is the same as Eq. (145ft . 

The stabilizer group of the graph code is given by 











B T Inn'-k' 



h' B 
-B T K 1 + W T -B T W + K 2 



(52) 



By Lemma [2] this matrix describes the same group as S c = ( M \ N ) because we have 









In'-k' 



o \ 




V 









Jn'-k' 



(n-l)n' 



-B T !„„/ 





nn'—k' 



If,/ B 

-B T K 1 + W T -B T W + K 2 



(53) 



□ 



Notice that from Eq. (|4ip . the adjacency matrix G c does not depend on the classical code of the 
outer code (C ou t), which indicates that Theorem Q] could also be true even if C ou t is nonlinear (in 
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this case Eq. f|42j) does no longer apply, but Eq. (|43p may still hold). As the stabilizer formalism 
can no longer be used to handle this case, we need an alternative proof technique which will be 
presented in the next section. This new technique is based on analyzing the encoding circuit of the 
concatenated code. It can easily be extended to more general cases, such as nonlinear outer codes, 
k > 1, and even to generalized concatenated quantum codes. 

V. GRAPH CONCATENATION BY GENERALIZED LOCAL COMPLEMENTATION 

In this section we prove our main result based on analyzing the encoding circuits of the con- 
catenated code. We start with an alternative proof for Theorem [1] in Sec. IV Al for the simple case 
that the inner code encodes only a single qupits, and the outer code is a graph code. This proof is 
based on the rule of "generalized local complementation." Then in Sec. IV B[ we show that the rule 
of "generalized local complementation" can be directly applied to the case that the outer code is a 
general CWS code, which is beyond the result of Theorem[T] In Sec. IV Cl we discuss the case where 
the inner code encodes more than one qupit (i.e., k > 1); we show that the rule of "generalized 
local complementation" given in Sec. IV Al also applies directly to this case, and hence completes 
the proof of the main result. 

A. Alterantive proof for Theorem 1 

Recall our main goal: suppose we have two graph codes Q out = (<5 ut>Cout) and Q m = (Qi n ,Cin) 
given by Eqs. (f33|) - (i37l) . where Q\ n encodes a single qupit. Let Q c = Q m C Qout denote the 
concatenation of the inner code Q m and the outer code Q ou t- We would like to show that Q c = 
(Qc,C c ), where Q c and C c are given in Eqs. ([41]) and (132]) . respectively. 

In Sec. [VAT] we first specify the encoding circuit of the concatenated code Q c , then we give the 
graphical interpretation of this circuit and define the encoding graph g^t{enc} ^ ^ e concatenated 
code Q c . Then in Sec. IV A~2l we define the rule of "generalized local complementation" on a graph; 
we show how the encoding circuit of the concatenated code Q c can be interpreted as generalized 
local complementation on the encoding graph, and how we can obtain the graph code Q^ c from the 
encoding graph g^t{enc}_ £ na rjy we that Q c and C c are exactly those given in Eqs. (I4ip and 
([42 p . thereby completing the proof. 

1. Encoding circuit and encoding graph for the concatenated code 

We have already discussed the encoding circuit of a concatenated code in Sec. Ill Bl Here we 
state it more formally. 

Procedure 3 (Encoding circuit for Q c with a graph outer code and an inner code encoding a 
single qupit) 

1. Apply the encoding circuit of Q ou t that encodes k' qupits into n' qupits which we call 
qi , . . . , q n i , as given by Procedure [1- 

2. Apply n' copies of the circuit that gives the graph state corresponding to Qi n . 

3. Apply W on all qupits q\, . . . , q n >. 

4- Apply the corresponding controlled-Z operators between these qupits and the graph states of 

Qin- 
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5. Apply H on q 1 , ...,?„/. 

6. Measure q±,. . . ,q n > in the computational basis. 
For an example, see the right circuit of FIG. El 

As discussed in Sec. Ill Bl Procedure [3] can be represented by a graph which is denoted by 
^Qc This graph is constructed as follows: Step 1 corresponds to the graph 5 7i Step 2 

corresponds to adding a copy of the graph Q m for each vertex qi of the graph <5 utj then we have a 
graph on n' + k' + nn' vertices; Steps 3, 4, 5 encode the n' qupits of the outer code into n' copies 
of the inner code, so we just add edges and labels according to controlled- Z gates that are applied 
between these n' qupits and the graph states of Q m . 

For an example of the encoding graph g^° ut i en6 ^ see ~piQ. [4j3. (The corresponding encoding 
circuit is given by the right circuit of FIG. [5j) 

2. Graph concatenation via Generalized Local Complementation 

We now give a graphical interpretation of Steps 3, 4, 5 given in Procedure El Notice that for 
1 < i < n! we apply to qupit qi, then the corresponding controlled-Z operations between qi 
and the z-th copy of the graph state Q- m , and finally H on qi. We show that each of these n' steps 
is equivalent to a generalized local complementation on the graph. 

Definition 1 (Generalized Local Complementation) Suppose F = (fij) is the adjacency matrix of 
a graph T , i is a vertex of J 7 , and fj is the i-th row of F. Also, let v be a vector whose coordinates 
are indexed by the vertices of J- such that v is zero on i and its neighbors, i.e., Vj = if j = i 
or fij 7^ 0. Then the generalized local complementation at (i, v) is the operation which sends F to 
F + v T fi + ff v. 

Notice that, since v is zero on the neighbors of i, for any j and k either (v T fj)jfc or (f/V)jfc is 
equal to zero. 

To get an idea on why we call this operation the generalized local complementation, let us 
consider the binary case. In this special case v corresponds to a subset of vertices (j belongs to 
this set iff Vj = 1). Then this operation is the same as to replace the bipartite graph induced on 
the neighbors of i and the vertices in v with its complement. (For an example, see FIG. [71) 

Theorem 2 (Encoding circuit interpreted as generalized local complementation) Consider a circuit 
which corresponds to a graph T with the adjacency matrix F . Let i be a vertex of T (or equivalently 
a qupit in the circuit), and let v be a vector which is zero on i and its neighbors. Suppose we change 
the circuit by applying on the i-th qupit, C z j (vj is the j-th coordinate o/vj on the qupits i and 
j, for any j, and then H on the i-th qupit. Then the resulting circuit is equivalent to the graph T 
after the generalized local complementation at (i,v). 

Proof: For simplicity assume i = 1, and let fi = (0 s), where fi is the first row of F. Also, 
let v = (0 v'), and J-' be the graph obtaining from T by removing its first vertex (and F' its 
adjacency matrix). Then the stabilizer group corresponding to the circuit is represented by 



s \ 
s T F'j 



(54) 
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Then based on the translation of the action of the Hadamard gate and controlled- Z gate on the 
stabilizer group (Lemmas [3] and H|), we can compute that stabilizer group after applying those gates 
as follows: 



s at 
s T F'j *U T I 



-1 s \ cl' v I 
F'j "~ L T / 



-1 

JT 



F' + s T v' 



H 



1 
-v' 



ST j 



s 
s T F' + s T v' 



(55) 
(56) 



Now to relate this stabilizer group to a graph, we change the set of generators by multiplying the 
above matrix by 



which gives 



1 

/ T 1, 



s 
s T F' + s T v' + v /T s 



(57) 



Hence, the adjacency matrix F of the graph is changed to F + v T fj + v. 
A direct corollary of Theorem 2 is the following: 



(58) 
□ 

C ut{enc} 



Corollary 1 Q c = (G C ,C C ), and the graph Q^ c can be obtained from the encoding graph £/g 
via Procedure^ 

Notice that Corollary Q] proves our main result in the case that the inner code encodes a single 
qupit and the outer code is linear. 

Also, note that the resulting graph of Corollary Q] is consistent with the one given by Theorem [T] 
since they both compute the same graph. In other words, the adjacency matrix of the graph Q^ c 
constructed via Corollary [T] is given by Theorem [TJ 

Theorem 3 The graph Q^ c given by Corollary[l\is equal to the graph given by Eqs. |^ i[ j and |^|j| ). 



Proof: Here we briefly describe a proof only for the binary case, and for the validity of Eq. (|4ip . 
This proof can simply be captured for the more general setting. 

Based on Procedure [lj the graph on which we apply the generalized local complementation 
operators has the following subgraphs: Q out with auxiliary vertices {1, ... ,n'}; and a copy of Q m 
with vertex set Vi for each auxiliary vertex 1 < % < n'. Then for each 1 < i < n' we apply the 
generalized local complementation on i with respect to Si Q Vi which is defined based on the 
classical inner code. 

Fact 1 Eq. describes the unique graph on the vertex set [J i Vi with the following structure: 

1. The induced subgraph on Vi, for every i, is isomorphic to Qi n . 

2. For every i ^ j, there is no edge between vertices in Vi and Vj \ Sj. 

3. For every i ^ j, there is an edge between vertices v G Si and w G Sj iff i and j are connected 
in Qout- 
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Clearly the graph with these properties is unique. Also, it is clear that Eq. (|1T|) represents this 
unique graph. 

Based on this fact, we will show that the graph resulting from Procedure Q] has the above 
structure. 

Fact 2 During Procedure [7] the changes on the subgraph induced on |J Vi happen only among ver- 
tices v G S{ and w G Sj for i ^ j. As a result, the final graph of Procedure^ satisfies properties 1 
and 2 of Fact El 

This is simply because in the generalized local complementations we never touch vertices in 
Vi \ S{. Furthermore, in each step, Si is disjoint from N(i) (the neighbors of vertex i), so there is 
no change in the subgraph induced on the vertex set Si. 

Fact 3 In Procedure El suppose we have applied the generalized local complementation on vertices 
1,2, ... ,1, for some 1 < I < n'. Then for any choice of Vi G Si, for 1 < i < I, the induced subgraph 
on vertices {vi, . . . , Vi} U {I + 1, . . . , n'} is isomorphic to Q ou t- 

This fact can be proved by a simple induction on I. 

Now we can prove the theorem. The resulting graph of Procedure Q] is a graph on the vertex set 
Uj Vi. According to Fact El this graph satisfies properties 1 and 2 of Fact [TJ Property 3 of Fact CD 
also holds based on Fact [3) Therefore, by the uniqueness of the graph described in Fact [H we are 
done. □ 

We illustrate the graph obtained by generalized local complementation for the code [[25, 1,9]] 
which can be obtained by self-concatenation of the code [[5,1,3]]. As a graph code, the code 
[[5,1,3]] can be described by a pentagon corresponding to the output nodes and a central input 
node that is connected to all output nodes. Using auxiliary nodes, the concatenated code [[25, 1, 9]] 
is show as the left graph in FIG. [H 



The outer code is given by the large pentagon with green/light dots and dashed lines. The five 
copies of the inner code correspond to the small pentagons with black dots and solid lines. The 
final graph is shown on the right of FIG. EE The five solid black pentagons remain, and any vertex 
in a small pentagon is connected with blue/dashed lines to any vertex of the neighboring pentagons 
as well as the central input node. 

The situation for the self-concatenation of Steane's code [[7,1,3]], which can be realized as a 
cube, is shown in FIG. [TOj 




FIG. 9: Self-concatenation of the code [[5, 1, 3]] yielding a code [[25, 1, 9]]. 
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FIG. 10: Self-concatenation of Steane's code [[7, 1, 3]] yielding a code [[49, 1, 9] 



B. A general outer code 

In this section we consider the case when the outer code is nonadditive. The advantage of 
Theorem [2] is that it directly applies to this case as well. 

Procedure 4 (Encoding circuit for Q c with a general outer code and an inner code encoding a 
single qupit) 

1. Apply the encoding circuit of Q ou t that encodes K' states into n' qupits which we call 
qi,...,q n >- 

2. Apply n' copies of the circuit that gives the graph state corresponding to Qi n . 

3. Apply W on all qupits q\, . . . , q n > . 

4- Apply the corresponding controlled-Z operators between these qupits and the graph states of 

Gin- 

5. Apply H on q u ...,q n >. 

6. Measure q±, . . . ,q n i in the computational basis. 
For an example, see the right circuit of FIG. [6l 

Notice that Theorem [2] deals with Steps 3, 4, 5 in Procedure [J which are exactly the same as 
Steps 3, 4, 5 as in Procedure [3j Therefore, whether C out is linear or not does not actually matter. 
Consequently, Corollary [JJ and thus the main result hold even for nonlinear outer codes. 



C. The case k > 1 

Theorem [2] can also be directly applied to the case when the inner code encodes more than one 
qupit. Again, to see this we only need to specify the encoding circuit of Q c . 

Procedure 5 (Encoding circuit for Q c with a general outer code and an inner code encoding k 
qupits ) 

1. Apply the encoding circuit of Q ou t that encodes K' states (or kk' qupits if C ou t is linear) into 
kn' qupits which we call qx,..., q^ n i . 
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2. Apply n' copies of the circuit that gives the graph state corresponding to Qi n . 

3. Apply H' on all qupits qi, . . . , qkn'- 

4- Apply the corresponding controlled-Z operators between these qupits and the graph states of 

Gin- 

5. Apply H on q u . . .,q kn >. 

6. Measure q\ , . . . , qtn' i n the computational basis. 

Note that Steps 3, 4, 5 remain the same as those given in ProcedureEJ Consequently, CorollarylU 
and hence our main result hold for the case of k > 1. 

For an example, the left graph of FIG. HUis the encoding graph g < ^ ut ^ enc ^ Q £ ^ e concatenated 
code Q c with a [[4, 2,2]]2 inner code, and a [[4, 2, 2]] 2 a outer code. Note that we decompose the 
outer code into two copies of a qubit code [[4, 2, 2]] 2 . Hence there are kn' = 2x4 = 8 auxiliary 
vertices (green/light vertices) in g^tienc} The corresponding encoding circuit is given by the 
middle circuit in FIG. \TT\ where "/" on each line indicates that there is a set of qubits, not just 
one. For instance, the line corresponding to \qo) represents the 4 input qubits (4 white vertices in 
the left graph of FIG. Qj] ), the line corresponding to \qi) represents the 8 auxiliary qubits, and 
the line corresponding to |<7 2 ) represents the 4 output qubits of a single inner code Q m . The graph 
Gc c of the concatenated code Q c can be obtained from the encoding graph g^° ut i enc J' by applying 
Corollary [TJ The result is shown as the right graph in FIG. [TT1 The blue/dashed lines are the 
edges obtained by generalized local complementation. 




FIG. 11: Graphs and encoding circuit for the concatenated [[16,4, 4]]2 code obtained by concatenating an 
inner code [[4, 2, 2]]2 with an outer code [[4, 2, 2]] 2 2. 



VI. GENERALIZED CONCATENATED CODES 

In this section, we discuss the application of our main result to the case of generalized con- 
catenated quantum codes (GCQCs). The construction of GCQCs has been recently introduced in 
Tol . 11). It resulted in many new QECCS, both stabilizer codes and nonadditive codes. 



A GCQC is derived from an inner quantum code = ((n, q\qi ■ ■ ■ q r , d%)) p , which is first 
partitioned into q\ mutually orthogonal subcodes Qj^iLi (0 < ii < <7i — 1), where each Qj^iLi is 
an ((n,q2 ■ ■ ■ q r ,d<2))p code. Then each SjjJiLi is partitioned into q2 mutually orthogonal subcodes 

(2) (2) 

^in{ui 2 } (° - i<2 ~ q2 ~ where each Sin{Hi 2 } has P arameters (( n ^3- • • qr,d 3 )) p , and so on. 
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f r _2) (r 1) 

Finally, each Q\ n i ili2 i 2 \ is partitioned into q r -\ mutually orthogonal subcodes Q- m i ili2 j 1 i = 
((n,q r ,d r )) p for < i r -% < q r -i — 1. Thus 

91-1 92-1 

2in = (J) 2in{ti}' ^in{n} = ^ 2 in {j li2 }, (59) 
«l=0 i 2 =0 

and d\ < d.2 < ■ ■ ■ < d r . In addition, we take as outer codes a collection of r quantum codes 
Qout' ■ ■ ■ ' Slut' where Q^ t is an ((n', K'-,d'-)) q . code over the Hilbert space %® n '. 
The generalized concatenated code Q 9C is a quantum code in the Hilbert space 

n ®nn ofdi _ 

mension K' = K'^Kl^ ■ ■ ■ K' r . The detailed construction of Q gc can be found in Here we only 
emphasize that the essence of the "generalization" , which is different from the usual concatenated 
quantum codes, is that the outer code is actually a product of r outer codes, and the inner code 
is nest-decomposed to specify how those product of outer codes are encoded into each inner code. 
Therefore, similar to a concatenated code Q c , a GCQC Q gc with a graph inner code 

sL 0) = (^ 0) ,4 0) ) (60) 

and r CWS outer codes 

s21 = (€,a (6i) 

naturally has an encoding graph, denoted by Gn nc \ and the corresponding encoding circuit is 
given by the following procedure. 

Procedure 6 (Encoding circuit for generalized concatenated code Q gc ) 



1. Apply the encoding circuits of Q^^t $iat encodes Kj states (or kjlog p qj qupits if C^t 
linear) into n'log p qj qupits which we call q±, ■ ■ . ,q n 'log p qj- 

2. Apply n' copies of the circuit that gives the graph state corresponding to Qi n . 

3. For each j = 1, . . . , r, apply on all qupits q\, . . . , q n 'iog p qj ■ 

4- Apply the corresponding controlled-Z operators between these qupits and the graph states of 
Qin- 

5. For each j = 1, . . . , r, apply H onq 1 ,...,q n , logp q . . 

6. For each j = 1, . . . ,r, measure qi, ■ ■ ■ , q n ' \ ogp q - in the computational basis. 

Notice that Steps 3, 4, 5 remain the same as those given in Procedure[3j Consequently, a similar 
result as in Corollary Q] holds for constructing GCQCs as well. 

Corollary 2 Q gc = (G gc ,C gc ), wh ere Q gc can be obtained for the encoding graph Qq ^ via Proce- 
dure U\ and C gc is the classical generalized concatenated code with inner code ( with correspond- 
ing decompos'i 
(j = I,-- - ,r) 



•is 



ing decomposition given by the decomposition of Q[„ , see [11] for details) and the outer codes C^ t 
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For an example, the left graph of FIG. [12] is the encoding graph Qq^ c of the GCQC Q gc 
with a [[4, 2, 2]]2 inner code code that is decomposed into two copies of a code [[4, 1,2]]2- There 
are two different outer codes [[4, 4, 1]]2 and [[4, 2, 2]]2- Note that there are 4 + 2 = 6 input vertices 
(white vertices) and 8 auxiliary vertices (green/light vertices). The corresponding encoding circuit 
is given by the middle circuit in FIG. \12\ where "/" on each line means that the line actually 
represents a set of qubits. For instance, the line corresponding to |goo) represents the 4 input 
qubits of the [[4,4, 1]]2 outer code, the line corresponding to \qoi) represents the 2 input qubits 
of the [[4, 2, 2]]2 outer code, the lines corresponding to \qio) and \qu) represents the 4 auxiliary 

qubits of the [[4,4, 1]]2 and the [[4, 2, 2]]2 outer codes, respectively, and the line corresponding to 

c 

| (72) represents the 4 output vertices in a single Qi n . To obtain the graph Q g 9 c c of the concatenated 

c' ' {enc} 1—1 

code Q gc from the encoding graph Gq u± apply Corollary [2J The result is shown as the right 
graph in FIG. [TTJ 




FIG. 12: Graphs and encoding circuit for the generalized concatenated [[16,6, 2]]2 code, derived from an 
inner code [[4, 2,2]] 2 and outer codes [[4, 2,2]] 2 and [[4,4, 1]] 2 . 



VII. CONCLUSION AND DISCUSSION 

In this paper we develop a systematic method for constructing concatenated quantum codes 
based on "graph concatenation" , where graphs representing the inner and outer codes are concate- 
nated via a simple graph operation called "generalized local complementation." The outer code is 
chosen from a large class of quantum codes, called CWS codes, which includes all the stabilizer 
codes as well as many good nonadditive codes. The inner code is chosen to be a stabilizer code. De- 
spite the restriction that the inner code must be a stabilizer code, our result applies to very general 
situations — both binary and nonbinary concatenated quantum codes, and their generalizations. 

Our results indicate that graphs indeed capture the "quantum part" of the QECCs. Once 
the graph part is taken care of, the construction of quantum code is reduced to a pure classical 
problem. This was essentially the idea of the CWS framework (i.e., the problem of constructing 
a CWS quantum code is reduced to the problem of finding a classical code with error patterns 
induced by a given graph). Here we have demonstrated that this idea extends to the construction 
of (generalized) concatenated quantum codes as well (i.e., to construct (generalized) concatenated 
quantum codes, given the rule of graph concatenation, one only needs to construct the (generalized) 
classical concatenated codes). We believe that our results shed light on the further understanding 
of the role that graphs play in the field of quantum error correction and other related areas in 
quantum information theory. 
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